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ABSTRACT 


There is a large gas density decrease from reactor inlet to exit in a nuclear rocket. 
Under certain reactor operating conditions, the momentum pressure drop can be a large 
portion of the overall pressure drop. As a result, deviations in flow rate or passage di- 
ameter can cause the reactor walls to overheat. A simplified analytical investigation has 
been made to find the sensitivity of the flow and wall temperature in heated passages as a 
function of Mach number for (1) a change in heat addition, and (2) nonuniform passage di- 
ameter. The Mach number range is from 0. 2 to 0. 9. 
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FLOW AND WALL-TEMPERATURE SENSITIVITY IN PARALLEL PASSAGES 
FOR LARGE INLET TO EXIT DENSITY RATIOS IN SUBSONIC FLOW 

by Meyer Reshotko 
Lewis Research Center 

SUMMARY 

There is a large gas density decrease from reactor inlet to exit in a nuclear rocket. 
Under certain reactor operating conditions, the momentum pressure drop can be a large 
portion of the overall pressure drop. As a result, deviations in flow rate or passage di- 
ameter can cause reactor walls to overheat. 

A simplified analytical investigation has been made to find the sensitivity of both flow 
and wall temperature in heated passages as a function of Mach number up to Mach 0. 9 for 
(1) a 6-percent increase in heat addition and (2) a 2-percent decrease in passage diam- 
eter. 

An analytical investigation has been conducted using reference conditions which are 
considered typical for nuclear- rocket- reactor applications. The reactor consists of 
fueled material with many parallel passages through which gas flows. One passage was 
studied with the assumption of a constant pressure drop across all the passages. The 
wall temperature is assumed constant. 

The results show that for a 6-percent increase in heat addition the sensitivities of 
flow and wall temperature with respect to Mach number are small. Neither sensitivity 
parameter increased by more than 10 percent. This same effect takes place for in- 
creases in heat addition of up to 12 percent. For a 2- percent decrease in passage diam- 
eter the sensitivities of flow and wall temperature with respect to Mach number are also 
small. Neither sensitivity parameter increased by more than 7 percent. Therefore, for 
large inlet to exit density ratios, the effect of momentum pressure drop on the flow and 
wall temperature sensitivities as one goes to increasing Mach numbers is small. How- 
ever, for a given Mach number of 0. 9 with the heat addition kept constant, a 3-percent 
decrease in passage diameter caused a 10-percent increase in surface temperature. 



INTRODUCTION 


A pressure drop occurs in any compressible flow system because of the forces nec- 
essary to overcome friction and to accelerate the fluid. These two components are gen- 
erally termed friction pressure drop and momentum pressure drop. Any decrease in 
density of the compressible fluid, such as heat addition, increases the momentum pres- 
sure drop if the flow is to be maintained. Heating the propellant in a nuclear -rocket 
reactor is one example of a situation where there is a large decrease in density of the 
fluid as it flows from inlet to outlet. Here, the pressure drop across the reactor core is 
established by the characteristics of many flow passages acting in parallel. Any de- 
crease in fluid density due to increased heat addition to one passage must be accom- 
plished by a decrease in flow rate through the passage. It has been shown (ref. 1), for 
example, that, if heat input is increased in one channel, laminar flow through parallel 
passages may be unstable. This could result in a continuous decrease in flow to the hot 
channel until failure occurs. Although the turbulent flow case is stable, some maldistri- 
bution will occur with nonuniform heating and can be expected to become more severe as 
the Mach number and momentum pressure drop increase. Maldistribution of flow can 
also occur with uniform heating but with nonuniform passage diameter. 

A simplified analytical investigation has been made to determine both flow and wall- 
temperature sensitivities in heated passages as a function of Mach number up to Mach 
0. 9 for a 

(1) 6-percent increase in heat addition 

(2) 2-percent decrease in passage diameter 

This report on flow and wall-temperature sensitivity is an analytical treatment using 
reference conditions which are considered typical for nuclear- rocket application. The 
wall temperature is assumed constant with respect to axial position. One flow passage 
of a multipassage configuration was studied with the assumption of constant pressure 
drop across all the passages. Simplified equations were developed to calculate the pas- 
sage operating conditions during steady state before and after a change was made in 
either heat addition or passage diameter. 


METHOD OF CALCULATION 

In analyzing the problem of flow and wall-temperature sensitivity the following reac- 
tor model was selected. The model consisted of fueled material with many parallel pas- 
sages through which the gas flows. The passages were all assumed to be uniform with 
respect to length, diameter, and heat addition except for one passage in which either di- 
ameter or heat addition were allowed to vary. This one passage was studied with the as- 


2 



sumption that the pressure drop was fixed by the reference operating conditions of the 
remaining parallel flow passages. 


Reference Conditions 

Reference conditions chosen to be typical of nuclear-rocket operation using hydrogen 
as the propellant- coolant were as follows: 


Surface temperature, T , °R; K 5000°: 2778 

Gas exit total temperature, T qu ^, °R; K 4500°; 2500 

Gas inlet total temperature, T\ n , °R; K 320°; 178 

Passage length, L, ft; m 3. 25; 0. 991 

Passage equivalent diameter, D, ft; mm 0. 0104; 3. 17 


The values given establish mass velocity for all reference conditions. In order to 
explore the effect of Mach number on flow and wall-temperature sensitivity with a fixed 
reference mass velocity, the outlet static pressure was allowed to vary to provide a 
series of reference conditions with outlet Mach numbers varying from 0. 2 to 0. 9. 

Flow and wall-temperature sensitivities were investigated by making incremental 
changes in the reference values of heat addition and passage diameter. Heat addition was 
increased by 6 percent and passage diameter decreased by 2 percent. This section de- 
velops simplified relations used to calculate the reference conditions and the conditions 
existing after incremental changes were made to these reference conditions. 

From the reference conditions, the mass velocity can be calculated. The total and 
static temperature, static pressure, heat-transfer coefficient, heat flux, and Mach num- 
ber profiles as a function of passage length for reference condition operation may then be 
calculated. 

Calculation of mass velocity . - The mass velocity is calculated by taking a heat bal- 
ance across an increment of passage length (see sketch (a)). 

The heat increase of the gas in the control volume is equal to the heat transferred to 
the gas at the wall. 
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( 1 ) 


GAC dT = h(T - T)ttD 61 
P s 


(Symbols are defined in appendix A. ) The surface temperature T is assumed to be 
constant along the passage length. 

The heat-transfer correlation used to calculate the heat-transfer coefficient was ob- 
tained from reference 2. This correlation which offers a good representation of experi- 
mental, turbulent heat-transfer data is expressed as follows: 


h = 0.021 - Re°‘ 8 Pr°" 4 
D 



-(0.2900+0. 00 19Z/D) 


( 2 ) 


All gas properties were evaluated at the bulk temperature from the data of reference 3. 
Substituting equation (2) into equation (1) yields 

^°* 2 p r °- 6 
0.084 

The reference mass velocity was obtained from equation (3) by trial and error. 
Equation (3) was numerically integrated from T = 320° R (178 K) at 1 = 0 to 
l = 3.2b feet (0. 991 m) using an assumed value of mass velocity. The integration was re- 
peated until a mass velocity was found which resulted in the reference outlet temperature 

of 4500° R (2500 K) at l = 3.2b feet (0. 991 m). The value of mass velocity found for the 

o 

reference conditions was 301 pounds mass per square foot per second (1466 kg/(m )(sec)). 
This value is independent of pressure level and therefore applies to all reference condi- 
tions regardless of outlet Mach number. 

Total temperature was recorded for the final integration and can be plotted as a func- 
tion of passage length. The local heat-transfer coefficient is calculated by evaluating 
equation (2), and the local heat flux by the equation 

q = h(T g - T) 

They too can be expressed as a function of passage length. 

Static- pressure d rop. - The static-pressure drop across the passage is found by in- 
tegrating a differential static- pressure drop from exit to entrance. The differential- 
pressure drop is 



(0. 2900+0. 0019Z/D) 


dT 


T s -T 


dZ 

D 


( 3 ) 
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(4) 


dp = -4f ^ 
D 



— dV 


The first term on the right represents friction pressure drop and the second term mo- 
mentum pressure drop. Numerical integration of the differential-pressure drop equation 
can be accomplished by writing all unknown quantities in terms of p, T, and l. (The 
relation between T and l is known from the numerical integration of eq. (3).) There- 
suiting equation for static-pressure drop is developed in appendix B. 


dp = 
dT 


0.092 
0. 084 



(o. 2900+0. oom/D) 



(-1 + E 1/2 ) 



1-2 (E -1 / 2 ) +—^~ (-1 + E 1/2 ) 

g c p 2 V - 1 


^(e- 1 / 2 )! 

g c P 


(5) 


The term E, is a dimensionless parameter defined as follows: 


E = 1 + 2 


y - G^R T 


3 c p 


A static-pressure profile is found by numerically integrating equation (5) from the 
exit to the entrance. 

At the exit, the static pressure is determined from the relation 


P 



t_ 

T 


( 6 ) 


where the values of G and T are the known reference conditions evaluated at the exit, 
and M is the desired exit Mach number. Corresponding to this Mach number, the isen- 
tropic static to total temperature ratio t/T at the exit is found from reference 4. Static 
pressure can then be plotted as a function of passage length. 

With both the static pressure and total temperature as a function of passage length 
known, a static temperature profile is calculated using the equation which converts total 
to static temperature. 
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( 7 ) 



/ \ 9 

1/2 

-1 + 

1 + 2[ y ' 1 J G R T 



V v ) s c p 2 



y - 1 G 2 R J. 


y g 2 
^ S C p 

The static pressure and static temperature as a function of passage length may be used to 
calculate a Mach number profile from the following relation: 

M = — + (8) 

P yyg c 

The total heat added to the gas is calculated using the equation 

« = GA Cp<T out - T ln ) (9) 

This completes the reference conditions for the typical nuclear-rocket model. 


Flow and Wall-Temperature Sensitivity 

Flow sensitivity for a change in heat addition (AG/G)/(AQ/Q) is defined as a frac- 
tional change in mass velocity resulting from a given fractional change in heat addition. 
Flow sensitivity for nonuniform passage diameter (AG/G)/(AD/D) is defined as a frac- 
tional change in mass velocity resulting from a given fractional change in passage diam- 
eter. Wall-temperature sensitivity for a change in heat addition (AT C /T)/(AQ/Q) is ex- 
pressed as a fractional change in wall temperature resulting from a given fractional 
change in heat addition. Wall-temperature sensitivity for nonuniform passage diameter 

(AT /T )/(AD/D) is expressed as a fractional change in wall temperature resulting from 
s s 

a given fractional change in passage diameter. 

The flow sensitivity and wall-temperature sensitivity was investigated for a 

(1) 6-percent increase in heat addition 

(2) 2-percent decrease in passage diameter 

The sensitivity terms are calculated at exit Mach numbers ranging from 0. 2 to 0. 9. One 
flow passage of a multipassage configuration was studied with the assumption of constant 
pressure drop across all the passages. A change in heat addition or passage diameter to 
one passage results in a change in flow rate through that passage, while the static- 
pressure drop and exit static pressure are constant. The static-pressure drop and the 
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exit static pressure are, however, different for each exit Mach number reference condi- 
tion. 

Equation (5) can be used to calculate the pressure drop between two points providing 
the surface temperature is known and is constant; the gas temperature between the two 
points can be expressed as a function of passage length. However, when used to deter- 
mine the flow and wall temperature sensitivity parameters, equation (5) will be one equa- 
tion of a set of simultaneous equations, the solution of which will necessitate a time- 
consuming iteration process. Therefore, an approximate formula for pressure drop, 
which will make future calculations easier, is used instead of equation (5). 

Equation (4) can be expressed in incremental form 


Ap = - 



( 10 ) 


The perfect gas law p = pRt and the definition of mass velocity pV = G may be used to 
express velocity as V = GR(t/p). By making the substitution 


V = 


V out + V in 


which is justified when the variation in velocity is approximately linear With passage 
length, equation (10) can be rewritten as follows: 

ap=-<A 

Sc 

For the large density ratios used herein 


|f L | out + ^in 


D 


\ p out 



( 11 ) 


-^2- >35 
Pout 

the term t. n /p. n can be neglected so that equation (10) takes the even simpler form. 


AP = “ 

(fi + l 

\ g 2 r 

*out 


V D 

’ s c 

Pout. 
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The friction factor f is assumed to be constant. 

Flow sensitivity . - Equation (12) is rearranged so that all the constant terms are on 
the left side and the terms affected by the change are on the right side. 


Ap 



. p out 
c 



= constant 


(13) 
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The product G t ^ is constant, therefore, t ^ varies inversely as G . A AG is ob- 
tained by subtracting the reference value from a new value of G. The new value of G 
was chosen such that the value of AG used was that value necessary to make a change in 
Q of approximately 6 percent. The new exit total temperature is then found using the 
relation: 


T out ^out + 


y 


9 

1 G^R out 


2 y g r 


2 

p out 


(14) 


The new values of G and T Qut are inserted in equation (9) to give a new value of Q. 
When the reference value of Q is subtracted from this new value, a AQ is obtained. 
This is the actual change in heat addition. The flow sensitivity term for a change in heat 
addition (AG/G)/(AQ/Q) is obtained in this manner. 

In a similar manner, the flow sensitivity due to nonuniform passage diameter may be 
found. This flow sensitivity parameter is calculated for a constant value of heat addition 
to the gas Q in addition to the condition of constant pressure drop across the flow pas- 
sage. 

Changing the equivalent diameter an amount AD allows the following three equations 
to be solved simultaneously for new values of G, t and T^: 


Ap 



^out 

R 


G 2 t 


out 


(13) 


T out *out + 


y - 1 G 2 R ^out 

27 g C D 2 

c p out 


(14) 
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(9) 


Q=G£D 2 C p (T out -T in ) 

A AG term is obtained by subtracting the reference value from the new value of G. 

Thus, the flow sensitivity for nonuniform passage diameter (AG/G)/(AD/D) is determined. 

Wall-tem perature sensitivity . - In determining the wall-temperature sensitivities 
due to a change in heat addition (AT /T )/(AQ/Q) and nonuniform passage diameter 
(AT /T )/(AD/D), a new surface temperature is found based on the new conditions. The 
new surface temperature is found by a trial and error technique using equation (3). In 
order to obtain a starting value for this technique, an approximate value of surface tem- 
perature was determined. From the heat balance of equation (1), the following equation 
is obtained: 


dT 

T s- T 


G C pD 


(15) 


Integrating this equation for the boundary conditions T = T in at 1 = 0 and T = T Qut at 
l = L the approximate surface temperature is 


T 


s 


T. 

m 


e 4(L/D)(h/GCp) T 


out 


1 - 


e 


4(L/D)(h/GCp) 


(16) 


The value of T ou j. has already been calculated in finding the flow sensitivity due to a 
change in heat addition. The heat- transfer coefficient h at the changed condition is 
found from the relation 



0 . 8 


where h ref and G re j are reference values. The term h rg ^ is found by evaluating 
equation (2) at the inlet and exit and dividing the sum by 2. With this as a starting value, 
repeated values of T g are used until equation (3) is balanced. After the new surface 
temperature is calculated, the difference between the new and reference surface temper- 
ature AT_ is determined. The surface-temperature sensitivity term (AT /T)/(AQ/Q) 
s s s 

is then calculated where AQ/Q is the same as in the calculation of flow sensitivity. 

A similar calculation is made for the wall-temperature sensitivity due to nonuniform 
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passage diameter (AT /T )/(AD/D). Surface temperature is determined in the same 
manner as in the previous sensitivity calculations except that the value used for T qu ^ is 
taken from the simultaneous solution of equations (9), (13), and (14) used in the calcula- 
tion of flow sensitivity due to a decrease in passage diameter. After the surface- 

temperature difference AT„ is determined, the surface-temperature sensitivity term 

s 

(AT /T )/(AD/D) is calculated where AD/D is the same as in the calculation of flow 
s s 

sensitivity. 


RESULTS AND DISCUSSION 

Methods for calculating reference conditions and conditions after changes, have been 
established in the section METHOD OF CALCULATIONS. Once the reference conditions 
have been established the sensitivity parameters may be determined for a change in heat 
addition and nonuniform passage diameter based on the developed simplified equations. 
Again, it should be emphasized that, during the changes, the static- pressure drop re- 
mains constant and the level of the wall temperature increases even though the wall tem- 
perature along the passage remains constant. 


Reference Conditions 

Using the iterative method described in the section METHOD OF CALCULATIONS 

along with the typical reference conditions for a nuclear rocket, the calculated reference 

hydrogen mass velocity is 301 pounds mass per square foot per second (1466 kg/ 

2 

(m )(sec)). The total temperature, heat- transfer coefficient, and heat-flux profiles are 
calculated and shown in figure 1. Exit static pressure and static-pressure drop are 
shown in table I as a function of exit Mach number. 

As an illustrative example for the case of M = 0. 9, the static pressure, Mach num- 
ber, and static temperature are plotted as a function of length and are shown in figure 2. 

For temperatures up to 2000° R (1111 K), the static temperature is almost the same 
as the total temperature. The static temperature reaches a peak of 4075° R (2264 K). 
This corresponds to a point 3. 05 feet (0. 930 m) from the entrance with a local Mach 
number of 0. 62. The static-temperature peak occurs at this Mach number because of the 
combination of heat addition and friction in the passage. For an ideal case without fric- 
tion, this peak occurs at Mach 0. 85 (ref. 5). 

The heat-transfer correlation used in equation (1) is taken from reference 2. It of- 
fers a good fit to experimental, turbulent heat-transfer data. Appendix C shows that dif- 
ferent heat-transfer correlations can markedly affect reference design conditions. 
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Position along length of passage, I, m 


Figure 1. - Total temperature, heat flux, and heat -transfer coefficient as function of passage length. Mass velocity, 301 pounds 
mass per square foot per second (1466 kg/(m 2 l(sec)l; surface temperature, 5000° R (2778 K); gas, hydrogen. 
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oT Sta H C temperature ' and Mach number as function of passage length for exit Mach number of 
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\tna is); gas, hydrogen. 


Mach number. 



TABLE I. - EXIT STATIC PRESSURE AND 


STATIC- PRESSURE DROP AS FUNCTION 
OF EXIT MACH NUMBER 


Mach 

number, 

M 

Exit static pressure, p Qu £ 

Static-pressure drop, Ap 

psi 

kN/m^ 

psi 

kN/m^ 

0.2 

2880 

19 857 

400 

2758 

.3 

1914 

13 200 

577 

3978 

.4 

1427 

9 839 

732 

5047 

. 5 

1134 

7 819 

867 

5978 

.6 

935 

6 447 

983 

6778 

.7 

792 

5 461 

1084 

7474 

.8 

681 

4 695 

1174 

8095 

.9 

600 

4 137 

1246 

8591 


Flow Sensitivity 

In figure 3 the flow sensitivity parameters (AG/G)/(AQ/Q) and (AG/G) /(AD/D) are 
plotted against Mach number. The change in flow sensitivity parameter due to an in- 
crease in heat addition as a function of Mach number is small, increasing from a value 
of -0. 81 at Mach 0. 2 to a value of -0. 91 at Mach 0. 9. These data are for an increase in 
heat addition of about 6 percent. The same effect was noted for several points calculated 
with an increase in heat addition of 12 percent. Thus, it appears that the percent 
change in mass velocity is nearly equal to the percent change in heat addition and this ef- 
fect occurs regardless of Mach number. The change in flow- sensitivity parameter due 
to a decrease in passage diameter as a function of Mach number is also small, increas- 
ing from a value of 2. 32 at Mach 0. 2 to a value of 2. 43 at Mach 0. 9. These data are for 
a decrease in passage diameter of 2 percent. 

The flow- sensitivity parameter does not vary greatly with Mach number because the 
expected increase in the ratio of momentum to friction pressure drop with increasing 
Mach number did not take place. The increase in this ratio did not take place due to the 
large density decrease of the gas from entrance to exit. 

For cases where the inlet density is not much greater than the outlet density, that is, 
Pi n /p ou t w 2, any decrease in the outlet density keeping the inlet density fixed, will cause 
an increase in the ratio of momentum to friction pressure drop. This is noted in the fol- 
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Figure 3. - Flow sensitivity parameters as function of Mach number. 


lowing equation for a given value of f(L/D): 


Ap = - 


cr 


Friction 


f L(_L + J_ l + 


Momentum 
1 1 


D 


V^out P in/ \Pout P i 


in/ 


(17) 


Equation (17) was obtained by making the substitution 


P 


= J_ 

Rt 


in equation (11). As the ratio Pi n /p ou t increases, the ratio of momentum to friction 
pressure drop increases until it reaches a point where any further increase in the density 
ratio will cause a very small increase in the pressure-drop ratio. For large changes in 
density which are typical for nuclear rockets, the term l/p out is greater than l/p in by 
at least a factor of 35. Thus, the term l/p in can be neglected and equation (17) takes 
on the following form : 


Ap (f-+l) (18) 

^c^out ' ^ ' 

This means that the ratio between momentum and friction pressure drop is a constant 
equal to the ratio l/f(L/D) irrespective of Mach number. Therefore, for large density 
changes there is no increase in momentum pressure drop relative to friction pressure 
drop for an increase in Mach number. 
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Wall -Temperature Sensitivity 


In figure 4, the wall-temperature sensitivity parameters (AT /T )/(AQ/Q) and 

b S 

(AT /T )/(AD/D) are plotted against Mach number. The change in wall temperature 
sensitivity parameter due to an increase in heat addition as a function of Mach number is 
small, increasing from a value of 1. 36 at a Mach number of 0. 2 to a value of 1. 46 at a 
Mach number of 0. 9. These data are for an increase in heat addition of 6 percent. The 
change in wall temperature sensitivity due to a decrease in passage diameter as a func- 
tion of Mach number is also small, increasing from a value of -3. 00 at Mach 0. 2 to a 
value of -3. 22 at Mach 0. 9. These data are for a decrease in passage diameter of 2 per- 
cent. 

Figure 5 is a plot of surface temperature T with respect to passage diameter D 

b 

for the reference Mach number of 0. 9 and with the heat addition Q kept constant. The 
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(a) Wall-temperature sensitivity due to change in heat addi- 
tion. 


tb) Wall -temperature sensitivity due to nonuniform passage 
diameter. 


Figure 4. - Wall-temperature sensitivity parameters as function of Mach number. 
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Figure 5. - Surface temperature as function of passage diam- 
eter for Mach 0.9. Constant heat addition; constant static 
pressure drop; gas, hydrogen. 
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passage surface reaches a temperature of 5500° R (3056 K), an increase of 500° R 
(278 K), for a decrease in diameter of 0.004 inch (0. 102 mm) or a 10-percent increase 
in temperature for a 3-percent decrease in equivalent diameter. 


Pressure Drop 


In the section METHOD OF CALCULATION, the pressure drop equation, equa- 
tion (12), is a good approximation to equation (5). 


Ap = - 



L 

D 



G 2 R t out 
g c ^out 


( 12 ) 


It uses a simple algebraic relation based on end conditions, to determine changed condi- 
tions for a constant pressure drop as a substitute for a complex equation which can only 
be solved by numerical means. The value of t t as found from equation (13) and based 
on equation (12), is correct to within a few degrees of the t .j. found by the exact method 
using equations (3) and (7). In a check made on the approximate solution, the flow sensi- 
tivity was found to be within 1 percent of that calculated using the exact solution. 

It should be noted that equation (12) was not intended to be used to evaluate pressure 
drop. If it were so used it would yield a maximum pressure drop error of 20 percent in 
the range covered by this report. 

In treating the flow sensitivity problem, constant static-pressure drop across the 
tube passages was used as an operating condition. One may ask why total- pressure drop 
was not used instead. The answer is that for this particular investigation there are large 
changes in density which allow the term l/p^ n in equation (17) to be neglected. Thus, 
whether constant static- pressure or constant total-pressure drop is assumed across the 
flow passages, the flow change will be influenced in the same manner. This can be 
shown as follows. Express equation (12) as a pressure ratio which gives 


^out 


G 2 R 

g c 




+ 1 


(19) 


Substituting the expression 
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into equation (19) yields a static- pressure ratio 


^=’' M out( f £ +1 ) + 1 < 2( » 

Pout V D ' 

For a given change the pressure ratio p in /p Qut is a constant due to the condition of 
constant- pres sure drop. The term y is constant, and, for the case of a change in heat 
addition, the term f(L/D) is constant too. For the case of nonuniform passage diameter, 
f(L/D) varies slightly due to the variation in D. But because its effect is negligible on 
equation (20), it is considered constant. Therefore, the term must be a constant. 

Now, consider the case of a change made at a constant total-pressure drop. Using 
the isentropic relations for the ratios of total to static pressure and equation (20) enables 
the ratio of inlet to exit total pressure to be expressed as 




in 




out 





y/y-i 


( 21 ) 


In this study the value of M in « 0. 06. Therefore, equation (21) can be approximated by 


Pin /P r 




r 


in ' p out 


in 


out 


rM out( 1 D + 1 ' + 1 



y/y-l 


( 22 ) 


Now the total-pressure drop is expressed as a ratio similar to the static- pressure ratio. 
By definition, & oxs x is constant in equation (22). The term y is constant, and, for 
the case of a change in heat addition, the term f(L/D) is constant too. For the case of 
nonuniform passage diameter, f(L/D) varies slightly, but, because its effect is negligible 
on equation (22), it is considered constant. Therefore, the term M qu ^ must be a con- 
stant. As a result, the term in the first bracket which is the static-pressure ratio 
Pi n /p ou t mu st also be constant. Furthermore, because M^ n «0.06, it may be as- 
sumed that p in = , because at that low a Mach number, the isentropic pressure ratio 

is almost 1. Therefore, p. is a constant and p . must also remain constant. Thus, 

’ *in *out 

the assumption of constant total- pressure drop is equivalent to the assumption of constant 
static- pressure drop for the case investigated in this report. 
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CONCLUSIONS 


A simplified parallel passage flow and wall-temperature sensitivity analysis was 
made for a typical nuclear-rocket reference condition. The conclusions resulting from 
changes in heat addition and passage diameter over a Mach number range of 0. 2 to 0. 9 
are as follows: 

1. For a change in heat addition with a constant static- pressure drop across the pas- 
sage, the change in the flow and wall-temperature sensitivity parameters (AG/G)/(AQ/Q) 

and (AT„/T )/(AQ/Q), respectively, are small with respect to Mach number. The flow 
s s 

sensitivity parameter increases from -0. 81 to -0. 91 and the wall-temperature sensitivity 
parameter increases from 1. 36 to 1. 46 over the Mach number range of 0. 2 to 0. 9. 

These data are for an increase in heat addition of 6 percent. This same effect takes 
place for increases in heat addition of up to 12 percent. 

2. For nonuniform passage diameter with a constant static-pressure drop across the 
passage and constant heat addition, the change in the flow and wall-temperature sensitvity 
parameters are small with respect to Mach number. The flow sensitivity parameter in- 
creases from 2. 32 to 2. 43, and the wall-temperature sensitivity parameter increases 
from -3. 00 to -3. 22 over the Mach number range of 0. 2 to 0. 9. The equivalent diameter 
is decreased by 2 percent. However, there is a large increase in surface temperature 
for a small decrease in passage diameter. At Mach 0. 9 a 3-percent decrease in diam- 
eter causes a 10-percent increase in surface temperature. 

3. For large inlet-outlet density ratios (>35) a formula was developed for approxi- 
mating pressure drop. This simplified pressure-drop formula is based on end conditions 
only. Using this simplified formula yielded a flow sensitivity within 1 percent of that cal- 
culated using the exact solution. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, March 11, 1968, 

120-27-04-56-22. 
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APPENDIX A 


SYMBOLS 


A 

2 2 

flow cross-sectional area, ft ; m 


2 2 

total pressure, lb force/ft ; N/m 

C p 

specific heat of gas at constant 

p 

2 2 

static pressure, lb force/ft ; N/m 


pressure, Btu/(lb mass)(°R); 

Q 

rate of heat transfer to gas, 


j/(kg)(K) 

Btu/sec; J/sec 

D 

passage diameter, ft; m (except 

q 

heat flux, Btu/ (sec) (ft 2 ) ; 


where otherwise noted) 


J/(sec)(m 2 ) 

f 

average friction factor 

R 

gas constant, (ft) (lb force)/ 

G 

mass velocity, lb mass/(see)(ft 2 ); 


(lb mass) (°R); (N)(m)/(kg)(K) 


kg/(m 2 )(sec) 

Re 

Reynolds number, DG/p 

g C 

gravitational constant, 

o 

T 

gas stagnation temperature, °R; K 


32.2 (lb mass)(ft)/(lb force)(sec ) 





T s 

surface temperature, R; K 

h 

heat-transfer coefficient, 

b 

gas static temperature, °R; K 


Btu/ (ft 2 ) (sec) (°R) ; J/(m 2 ) (sec) (K) 

t 

k 

thermal conductivity of gas, 

V 

gas velocity, ft/sec; m/sec 


Btu/ (ft) (sec) (°R) ; J/(m)(sec)(K) 

r 

ratio of specific heats of gas 

L 

length of passage, ft; m 

M 

absolute viscosity of gas, 

l 

position along length of passage, 


lb mass/(sec)(ft); (N)(sec)/m 2 


ft; m 

P 

gas density, lb mass/ft 2 ; kg/m 2 

M 

Mach number 

Subscripts: 

Nu 

Nusselt number, hD/k 

b 

bulk 

P 

diametral perimeter of passage, 

in 

inlet 


ft; m 

out 

exit 

Pr 

Prandtl number, Cp/i/k 

ref 

value at reference point 



s 

surface 
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APPENDIX B 


DERIVATION OF PRESSURE DROP DERIVATIVE WITH RESPECT TO TOTAL 

GAS TEMPERATURES 


Pressure drop can be expressed as the sum of the friction and momentum pressure 
drops 


dp = -4f ^ dV (Bl) 

D 2 §c g c 

The definition of mass velocity G = pV is substituted into equation (Bl) to form the equa- 
tion 


dp = 


lf & GV 
D 2g c 


— dV 
g c 


Using the definition of mass velocity G, and the perfect gas law 


P 


_ _P_ 
Rt 


results in an expression for velocity 


V = GR- 
P 


The equation for total temperature 


(B2) 


(B3) 


T = t + - 
2 


1 (y ~ A g 2 r t 2 


g 2 
& c p 


may be rearranged to express static temperature as 
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/ \ 2 1 

1/2 

-1 4- 

1 + 2[ y “ 1] G 2 R T 

V y i Sc P 2 



y - 1 G 2 R 1 


y g 2 
r & c p 


Gas velocity can now be expressed in terms of mass velocity, total gas temperature, and 
static pressure. 


GR<-1 + 


1 + 2 (^ T 
V 7 ; 


11/2 


V = 


c p 2 J 




(B4) 


y - 1 G R 1 

y g c P 


The velocity differential with respect to the variables T and p is 


r 


dV = GR/ e 

y - 1 G 2 R 

y g c 


d 


-1 + 




f 2 V / 2 

•1 + i + 2 rjLl&il.) 


y 


g 2 
& c p 


V 


-I y 


dp 

J. g 2 r[ 


(B5) 


The first term in equation (B5) can be simplified to 
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2 ' 1/2 

.i + |i +2 r^i^5X 

^ g C p 2 


2 / 2 \ _ 

= y - 1 G_R L + 2 y - 1 G^R _T \ dT 


y s r 


y g 2 

f & c p 


2 / 2 \“ 

_ 2 yjii g^r [i + 2>l^ g^r JL) T dp 


y g„ 


y g 2 
r & C p ' 


Substituting it back into equation (B5) gives a more simplified velocity equation. 


r 


■ 1/2 

dV = GR^| 1 + 2 31ZJ: G!RX] dT- 2 [l + 2>^i^ T 

p 


rV2 


y g 2 
r & c p ' 


v. 


y g c p 2 j 


dp 


'N 


f 2 X 1 / 2 

-l + [l + 2 IliG!Rl 

77 g c p 2 


dp 


V (B6) 


y-1 G R 

y g n 


Solve equation (1) for dZ 


GACp 

hP 


dT 

Tg _ T 


(B7) 
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In equation (B8) the following term can be simplified by making the substitutions 
A = (tt/ 4)D 2 and P = jtD 


fGAC , GC 
P _ £ P 

hPD 4 h 


(B9) 
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The Stanton number for this particular correlation (ref. 2) can be expressed as 


h _ Nu 
GCp RePr 


/ T \ - (0. 2 900+0. 00 19Z/D) 
0. 021 Re 0, 8 Pr°- 4 (— j 

RePr 


A turbulent friction factor 


f = 


0. 046 


Re 


0.2 


taken from reference 6 changes equation (B9) to 


fGAC p 0. 046 
hpD "" 4Re°* 2 


RePr 


/ T \ -(0. 2900+0. 0019Z/D) 
0. 021 Re 0, 8 Pr°- 4 (— j 


/„ \ (0. 2900+0. 0019Z/D) 
0- 046 P r°. 6j _s j 

0. 084 It j 


Inserting this into equation (B8) gives a formula for differential-pressure drop 


/„ \ (0. 2900+0. 0019Z/D) , 

0.092 ^0. 6/M G z R 

0.084 


dp = - 


Sr 


+ 

t r 
5 c p 


\l/2 


V S, r ?J 


dT 


T s ” T 


2 / 2 V 1 / 2 

G R ( i 1 2 y - 1 G R T ' 


y - 1 G z R 

r 


S„ 


y g 2 

y & c p > 



g 2 r 

-1 + (l + 2 - 

0 . 

( 2 Y 1/2 „ gc 

V 

G Z R 

i . o y - l g z r t \ _ dp 


T 

& c 

1 + Ct 11 

V V g c p 2 / P 2 

y - 


a/a 


y e c p 2 y 


dp 


(B10) 
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Rearranging terms in equation (BIO) yields the derivative of pressure with respect to 
total gas temperature 


0.092 pr 0.6(^s 
0. 084 


4R = 

dT 


(0.2900+0.00191/D) 


H 


g 2 r _t 

g c p 2 


1/2 


p G R 


T - T 




g 2 r _t 

g 2 
s c p 


- 1/2 


1-2 


G 2 R 


/ \ 2 

-1/2 

r 

/ \ 2 

1/2 

1 Jy - 1\ G Z R T 

JL +_>!_; 

-i+ 

i , 2 (r - l\ g R t 


V v ) Sc p 2 _ 

p 2 1" 1 


UK p 2 



(Bll) 


Equation (Bll) can be further simplified by the substitution 


1 + 2 


'y - 1) G~R T _ 


= E 


’c p 


The derivative of static pressure with respect to total gas temperature is 


dp = 

dT 


0.092 pr 0. 6/ 1 s 
0.084 


' T \(0. 2900+0. 0019Z/D) 

C* * 

T 


/ — 1— ) (-1 + E 1 ^ 2 ) P - G?R (e-1/2) 1 

\y - v T s - T g r p 


i- 2 g!R(e- 1 / 2 ) t. + _l_ ( -i + e 1 / 2 ) 


r - i 


The relation between T and l is known from the numerical integration of equation (3) 
of the text. 
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APPENDIX C 


DIFFERENCES IN HEAT-TRANSFER CORRELATIONS 


Heat-transfer calculations can be no better than the correlation equation which is 
used. One difficulty is that of choosing a heat-transfer correlation that best represents 
the conditions under investigation. Different heat-transfer correlations can markedly af- 
fect reference design conditions. As an example, a correlation based on properties eval- 
uated at bulk temperature with an l/D effect 


Nu b = 


h b D 




-(0.2900+0. 0019Z/D) 


(Cl) 


is compared with a correlation based on properties evaluated at surface temperature 


OR 04 

Nu = = 0. 021 Re~- e Pr^ * 

b S S 



■ 0.8 


(C2) 


For typical reactor operating conditions, keeping tube length, equivalent diameter, tem- 
perature and pressure levels the same, a bulk correlation will give a flow rate almost 
twice as large as that of a surface correlation. 

Studying the two correlations as they appear in the following two heat balances ex- 
plains this difference. The correlation which evaluates properties at the bulk tempera- 
ture and has a variable exponent with length in the T /T term is 

s 


G 


0.2 



(0. 2900+0. 00191/D) 


dT 


0.084 


Tg - T 


dZ 

D 


(C3) 


The correlation which evaluates properties at the surface temperature and has a 

constant exponent in the T /T term is 

s 
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(C4) 


G 0 ' 2 




0.084 



dT 

t s - T 


= dz 

D 


The left-hand differential terms in equations (C3) and (C4) are equal to the same in- 
crement dZ/D. The bracketed term in equation (C4) is larger than that of equation (C3). 
This occurs because the exponent in equation (C4) is a constant equal to 0. 8, while the ex- 
ponent in equation (C3) is a variable whose lowest value is 0. 29 and only exceeds 0. 8 near 
the end of the tube. Because both left-hand terms are equal to the same increment dZ/D, 

the mass flow rate in the bulk correlation (eq. (C3)) is greater than that in the surface 

0 2 

correlation (eq. (C4)). It must be remembered that small changes in (GD/p) ' cause 
large changes in GD/ju. 

The following table compares the results of the two: 


Correlation 

(GD/fi) 0 ' 2 

GD/fx 

Mass velocity, G 




lb mass/(ft 2 )(sec) 

kg/ (m 2 ) (sec) 

Bulk 

10.4 

122 000 

301 

1466 

Surface 

9. 1 

62 400 

171 

833 


0 2 

A bulk correlation will have a 14- percent increase in (GD/ju) ' which causes a 
76-percent increase in flow rate. 

The bulk correlation with the variable exponent, modified Dalle Donne correlation, 
was used because it offers the best agreement with experimental heat-transfer data (see 
ref. 2). 
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